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Building on previous works of the mid 1960's, we construct an integral equation for forward 
elastic scattering (t = 0) at arbitrary virtuality Q 2 and large s = W 2 . This equation sums the 
ladder production of massless intermediate bosons to all orders, and the solution exhibits Regge 
behavior. The equation is used to study scattering in a simple x x4> scalar theory, where it is solved 
appoximately and applied to the study of DIS at small x. We find that the model can naturally 
describe the quark distribution in both the large x region and the small x region dominated by 
Reggeon exchange. 



I. INTRODUCTION 

This paper discusses the ingredients of a theory of deep 
inelastic scattering (DIS) that holds for all values of the 
Bjorken scaling variable x — Q 2 /2Mbv (where Mb is the 
mass of the target bound state). The Q 2 of the virtual 
photon is fixed at a large value, so the Q 2 evolution of the 
structure functions will not be discussed. The individual 
ingredients are well known; our emphasis here is on com- 
bining these into a unified picture and in explaining how 
the various pieces fit together. 

At large Q 2 , the QCD coupling a s is small, and the 
strong interaction dynamics can be treated perturba- 
tively (i.e. to lowest order) except in regions where contri- 
bution from higher order loops are enhanced, so that con- 
tributions of all orders must be summed. It is well known 
that the existence of the bound state itself (the target qQ 
system in the simple example discussed in this paper) de- 
pends on the enhancement of ladder type loops, so the 
description of the bound state vertex function must in- 
volve nonperturbative physics [l8j. Nonperturbative con- 
tributions are also known to be important in the region 
of small x [1 H H II . In the 1960's it was shown that 
Regge behavior is generated by infinite sums of ladder di- 
agrams in the i-channel and explicit equations have been 
written down for these ladder- type amplitudes 0, ||| and 
solved analytically and numerically [6]. Furthermore, it 
was shown in the case of DIS that the finite scaling limit 
in x of structure functions with Regge behaviour is pos- 
sible to obtain if the Regge residues are properly tuned 
[7]. In what follows we derive explicitly the Regge be- 
haviour of structure functions in our model and show that 
it agrees with general expectations [1-4] and the conjec- 
ture made in Ref. 0. 

These basic ideas, and the specific model discussed in 
this paper, are illustrated in Fig. ^ The two diagrams 
show contributions to DIS from a bound state of a heavy 
anti-quark (or diquark, if applied to baryons) Q of mass 



M (assumed to have no charge, for simplicity) and a light 
quark q of mass m, interacting through the exchange of 
scalar gluons ((/> particles with mass y,; confinement is 
neglected in this discussion). To simulate some aspects 
of the behavior of QCD at high Q 2 , the "quark- gluon" 
coupling, g, will be assumed to be small. Even for small 
couplings bound states exist (as in atomic physics, for ex- 
ample) ; they emerge from the infinite sum of ladder- like 
diagrams in the bound state s-channel @- The second 
diagram, Fig.^b), describes contributions from the ex- 
pansion of higher qQ(f> n Fock states, and would normally 
be suppressed by the smallness of the coupling constant 
g. These diagrams are indeed small, except in the region 
of large Q 2 and small x, where each of the loops in the 
t-channel ladder sum is large, and contributions of this 
type must be summed to all orders. This nonperturba- 
tive contribution was discussed in Ref. 6] , and shown to 
exhibt Regge-like behavior. 

In accordance with Regge theory [1-4,7], at small x the 
quark distributions F\ and F 2 should scale as 



Wi{Q z ,x) -> F 1 (x) ~ Ai ( - 

1 x 



vW 2 (Q 2 ,x) -> F 2 (x) ~ A 2 
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where £(0) is the intercept of the corresponding Regge 
trajectory. The valence quark distributions f(x) = 
F2(x)/x are known to behave at small x as x~°- 5 while 
for sea quarks this behavior is a; -1 , corresponding to the 
Regge trajectories of the p meson and pomeron respec- 
tively Q . However the direct matching of the Regge pole 
contribution to the scaling behavior at larger x requires a 
specific tuning of the coefficients A.; , which is not known 
on general grounds [2] . The relevant behaviour of these 
coefficients was suggested in Ref. based on the avail- 
able solutions of scalar ladder equations with a special 
relation between the masses or upper bounds for the 
massive case US* . It would be instructive to use one solv- 
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FIG. 1: The two contributions to DIS studied in this paper, 
(a) Pole contribution, important at medium to large x. (b) 
Regge contribution arising from the contributions of higher 
Fock components, which dominates at small x. The expected 
qualitative behavior of each contribution is shown below each 
diagram. 



able model with Regge behavior to understand how the 
small x behavior explicitly matches with the behavior at 
x around one. 

This is the primary purpose of the present paper: to 
show how these two processes can be added together to 
produce a description of DIS valid for all x. We will carry 
out this discussion using the covariant spectator equation 
in the ladder approximation [Tl|. 

Our treatment of the diagram ^b) is based on the in- 
tegral equation for the imaginary part of the amplitude 
for the scattering of two scalar particles of mass m in 
the scalar % 2 </> theory, derived independently bysimilar, 
but slightly different methods, in Refs. p| and Q. From 
this equation, which essentially sums up ladder-type di- 
agrams, the intercept of a Regge trajectory was found 
analytically. It appears that this integral equation is a 
very useful tool for studying structure functions, as will 
be seen in what follows. 

This paper is divided into five sections. In the next 
section we derive the integral equation for forward "fq 
elastic scattering, allowing the external photons to have 
large Q 2 , which is relevant to applications to DIS, and 
the quark to have a virtual mass p 2 ^ m 2 , needed for 
the application to bound quarks. In Sec. Ill we solve the 
equation, and in Sec. IV we add the photon spin (which 
was neglected in Sec. Ill for simplicity), and explicily 
evaluate the two diagrams shown in Fig. ^ Our conclu- 
sions are presented in Sec. V and the Appendices include 
some further details. 

This work can be easily extended to the description of 
non-forward scattering, and can be used to model Gen- 
eralized Parton Distributions (gpd's) [HEIH 

which 

is one of the motivations of the present work, but these 
applications are saved for later papers. 



II. INTEGRAL EQUATION FOR 
PHOTON-QUARK SCATTERING 

The integral equation to be derived is illustrated in 
Fig. |21 We have specialized the kinematics to forward 
scattering, appropriate to applications to DIS, and as- 
sume here that the photon has zero spin (this will be 
corrected in Sec. IV). (A more general equation applica- 
ble at momentum transfers t ^ can be derived using the 
methods of Refs. 0,0, and will be discussed elsewhere.) 

The integral equation can be easily obtained from 
study of the box diagram, Fig. |3] The imaginary part 
of this amplitude is the total forward cross section for 
the process 7 + q — > q + g. Hence 



ImT 2 = 



2 2 

e g 



eV 
16tt 2 



3 p {2ir) 4 5 4 (k+p' -q-p) 



(2Tr)HLu k E p , 
d 3 kS(E p / + lou 



Po 



{p-k) 2 } 2 
- <7o) 



uJkEpi [m 2 - (p - k) 2 



(3) 



where e is the charge of the quark, and u>k = y /i 2 + k 2 
and E p i = \/m 2 +p' 2 . Since q 2 = —Q 2 < 0, it is conve- 
nient to evaluate this integral in the Breit frame, where 



q 
p 



(0,0, 0,Q) 
(po,0,0,^) 



(4) 



[The external quark is off-shell, and its energy is fixed by 
other variables.] Then the integrand is independent of 
the azymuthal angle, and the delta function may be used 
to evaluate the integral over k 2 ^, giving 



lm.T 2 (s,x,Q 2 



e 2 g 2 I ' dk 



57T 



ik* Po [m 2 - (p-k) 2 } 2 
where we have introduced the variables 

s = (p + qf = P 2 + Q 2 ( - - 1 J > (m + /i) 2 



(5) 



Q 2 

2p-q 



2p z 



(6) 





(a) 



FIG. 2: Diagrammatic representation of Eq. [\'2'M . The solid 
lines are "quarks" of mass m, the horizontal dashed line is a 
"gluon" of mass fi — * 0, and the vertical dashed line reminds 
us that all intermediate particles are on shell. Fig.(b) is the 
pole term. 
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Note that, because the external quark is off-shell, x is 
not restricted, and 



(7) 



The limits on the k z integral will be discussed below. 

It is convenient to write the internal integration vari- 
ables in terms of s' and x' where 



s' = (p-k + qf = m 2 = (p - k) 2 + Q 2 ( - t - 1 

1 X 



Q 2 



Q 



2(p -k)-q 2k z + § ' 

This substitution gives the result 

lmT 2 (s,x, Q 2 ) 

_ e 2 g 2 Q f x + dx' {s + rj 2 )- 1 ' 2 



(8) 



16tt J x x' 2 [m 2 + Q 2 - 1) - s'] 



'12 



with 



n = Q i 



2x 



(9) 



(10) 



The limits on the x' (or & z ) integration follow from the 
conservation of energy and the requirement that k\ > 0. 
The energy relation is 

y/fj, 2 + k\ + k 2 + ^s' + fc 2 + (t? - k z f = p Q . (11) 
This is a quadradic equation in k z , with roots at 



/.■: l (/a.) = //.Jd • /•-' 



(12) 



where 



.4 



s - s' + /i 2 
2~s 




np-k 




FIG. 3: Diagrammatic representation of box diagram, Eq. 10. 
The lines are as described in Fig. [2] 

The restrictions on s and s' insure that both A and A 
are real and positive. Introducing k f = kf(0), we see 
that as long as 



(14) 



there always exists a value of k\ > that satisfies Ijllf) 
and the integral is nonzero. This determines the limits. 
The corresponding limits on x' are 



x± 



1 + (2xkf)/Q 



1 + (2x-1)At^V(s+F)A 



(15) 



It is now a starightforward matter to turn the result 
(PJ) for the box diagram into an integral equation for the 
ladder sum. First observe that the single quark pole term 



ImTr^Q 2 ) = e 2 J (2nfS i (p'-p-q) 
= 2ire 2 6(s - m 2 ) , (16) 



A = [s — (vV + ji) 2 ] [s — (y/P — fi) 2 ] . (13) With this notation the box diagram can be rewritten 



9s 9 2 Q f^-ri (& [ x+ dx' ImTi(s',x',Q 2 ) , N 

ImT 2 ( S ,x,Q 2 ) = |-^ / -p=/ -To T-o ; \V I7T^ 



32tt 2 7 m2 yJT^ 2 y*_ ^' 2 K + Q 2 - 1) - a'] 



where the upper limit on the s' integral follows from the general requirement that *Js > \/P The integral equation 
we seek is the immediate generalization of (II 71) : 

t t ( rfl\ ww rfl\ m 9 2 Q f (V ~ S ^ )2 ds ' [ V+dx ' ImT^s^x^Q 2 ) 
It is instructive to cast this equation into an alternative form convenient for solution. To this end, first introduce 
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virtual mass parameters for the external and internal quarks: 



1 



- 1 



( P -k) 2 = Q 2 



For simplicity, we also set fi = 0. It is then a straightforward matter to transform Eq. i|18|) into 



(19) 



ImTo^s, k, Q 2 ) = lmTx(s, «, Q 2 ) + 



16tt 2 



ds' 



[to 2 + k'} 2 



(20) 



with 



2.s 



Q 2 + s - k± vA 



A = s 2 + k 2 + Q 4 - 2kQ 2 + 2sk + 2sQ 2 



(21) 



Finally, introduce new variables 



XB 



Q 2 



s' + Q 2 



s + Q 2 



In terms of these, Eq. 1|20(1 becomes finally 



(22) 



M{xb, k, Q 2 ) - M po \ c (xB, k, Q 2 ) 



2n2 /-i /•«+ dK > M(x B /a, K ',Q 2 ) 
da ' 



g 2 Q 



16n 2 xb 



K _ VS^ [m 2 + K >] 



(23) 



with 



A. Approximate equation 



On 



xb 1 



Q 2 

1 - a 



K± = K 



An = 



2(1 -x B ) 



xb 



- K ± VA 



k H — ■ — I - 4kQ z , 



with ImTo^s, k, Q 2 ) = M(xb, k, Q 2 ), and 

/ ^~*2\ 27TG X d c , 

M polc (x B ,K,Q ) = — — — d \ x B 



Q 2 



Q 2 



m 2 + Q 2 



(24) 



(25) 



Equation (|18fl is exact, and Eq. I|23|l is exact in the fj, = 
limit. 



The form of the r.h.s. of 123|) suggests that it will be 
large only if xb is small, and hence this term may be ap- 
proximated by taking xb small. This will be assumed in 
what follows, and confirmed later, after the approximate 
solution has been obtained. As -> 0, the following 
approximations hold 



7 A 

K+ 



xb 

Km; 



K + 0(X B ) 



+ 0(x B ) 
a k + 0(xb) 



l 



X B 



Q ( — + 1 ) (1 - a) 



(26) 



III. THE REGGE LIMIT 

We now apply Eq. H'-'il) to the case of DIS from an 
off shell quark (neglecting the photon spin until the next 
section). If the strong coupling g is small, under what 
circumstances (if any) will the r.h.s. be larger than the 
pole term, requiring a nonperturbative solution to the 
equation? 



This gives the following equation, which should hold ap- 
proximately for all x b and Q 2 >j> m 2 

M(x B , k,Q 2 ) - M po i c (x B ,K,Q 2 ) 



9 



16tt 2 



i 

da 



dK , M(xB/a,<Q 2 ) (2?) 



\ni + k 



l\2 



We emphasize that this equation differs from the exact 
Eq. (|23[l only in the form of r.h.s., which by assumption 
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is large only at small xb where the approximate form 
used in (|2*?|) is accurate. 



this anzatz into i|27|) . and neglecting the inhomogeneous 
pole term which is small at small xb, gives an equation 
for R(xb 7 k, Q 2 ) 



B. Approximate solution for small xb 

Following Ref. , assume that M exhibits a Regge-like 
behavior at small xb 



M(x b ,k,Q 2 ) 



where R is finite as xb — > 0. The power £ will be deter- 
mined self consistently from the equation. Substituting 







1 small xb 


{xb) 



R(x B , k, Q 2 ) 

1 (lag* /^V R ^ B l a ^f\ (29) 



16tt 2 



X B J OtK 



\wr + K 



where, by assumption, R(0, k, Q 2 ) is finite. This quantity 
satisfies the equation 



i?(0,K,< 



A to / da a 

lo J<u; 



dn 



, R(0,k',Q 2 ) 



dn! 
o V K 



- 1 ! R(Q,k',Q 2 ) 
| [to 2 + n'} 2 



(30) 



where A = g 2 /(167r 2 m 2 ) is dimensionless. Note that 
R(0, k, Q 2 ) is independent of Q 2 . The boundary con- 
ditions on R are R — > as n — * oo, and i?(0, 0,Q 2 ) = 
constant. Introducing the dimensionless variables y = 
k/tti 2 and y' — k' /to 2 , and differentiating the equation 
by y twice, leads to a differential equation for r(y) = 
R(0,k,Q 2 ) 

yr " + (£ + 2)r'+ * r = 0. (31) 

The general solution of this equation is a hypergeometric 
function Gj| 



r{y) oc F [ a,b,c, 



y 



l + y 



(32) 



with 



1 1 i— 

- 2 -- 2 VT 

1 + 2. 



4A 



4A 



(33) 

The value of i is fixed by the y -— > oo boundary con- 
dition, which requires that the hypergeometric function 
approach zero as z = y/(l +y) — > 1. First, the hypergeo- 
metric series will not converge at z = 1 unless c—a—b > 0, 
and this requires i + 1 > 0. Then, using the properties 
of the hypergeometric function |1RI |. we have 

T(c)T(c - a-b) 



F(a,b,c,l) = 



T(c-a)T(c-b) 
T{£ + 2)T(£+ 1) 



T(£- 



; VT+4A)r(£- 



4A) 



(34) 



This can be zero if and only if the arguments of one of 
the r functions in the denominator is —n, where n is 
zero or a positive integer. Since I + 1 > 0, it must be the 
argument of the second T function that is —n, and hence 
the only possible values of £ are 

1 



£ = £ n = - (^\/IT4A-3-2n) n = (0, 1, 2, • • •) . (35) 



The requirement £ + 1 > yields 

. (2£ + 3 + 2n) 2 - 1 n 
A = : > 0, 



(36) 



so there are solutions even if A is very small, as antici- 
pated. The solutions for n — were previously obtained 
in Ref. @. 

For large A there will be more than one value of I 
that satisfies the condition l|35|) . The number of values 
depends on the size of A. Specifically, 



£ = £ a = \ (VI + 4A - 3) if 2>A>0 
£ = t or £ x if 6 > A > 2 



(37) 



t = {t Q ,h,ti} 



if 12 > A > 6 



For A > 2 the general small x b solution will be a sum of 
contributions from the different trajectories 



M(x b ,k, Q 2 )\ 



;mall 



n=0 v 



(38) 
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The solution ro(y) for the "leading" trajectory, n — 0, is 

For small A (i.e. A < 2) there is only one trajectory, and 
we obtain the unique approximation 

/ 1 V°f m 2 V a+1 
M (* B ,^) = co(_) , (40) 

with c an undetermined coefficient. 

Some additional arguments which demonstrate that we 
have found all of the solutions are given in Appendix A 

At this point we are left with several unanswered ques- 
tions: (i) How is the coefficient Co to be fixed? (ii) Only 
the leading dependence at small x B is given by Mq. How 
large are the nonleading terms at small xb and what is 



their analytic form? (iii) How is the pole term to be in- 
cluded in the solution and how is the large x b dependence 
described approximately? 

All of these questions will be answered in the following 
subsection, ft turns out that substitution of the ansatz 

M(x B , k, Q 2 ) = M pole (x B ,K, Q 2 ) + M (x b ,k, Q 2 ) , (41) 

into Eq. H27|) will show how to fix the coefficient Cq so 
that (|41|l will give a satisfactory solution to 1|27|) for all 
x B . 



C. Approximate solution for all xb 

Substituting Mq into the r.h.s. of the asymptotic Eq. 
(|2Tjl gives 



r.h.s. 



i / \ to 
a 



coA m da [ — 



c A 



da 

XB 



( 



a 
\xb 



K [to 2 + K'] e °+ 3 

lo+2 



(4 + 2) 

The first term is the leading term, and doing the a integration gives 



lo+2 



(42) 



r.h.s. 



c A 



1 



(4 + 2)(4 + l) a4° \m 2 + n 



M (x b ,k, Q 2 ) - c x B 



xbk 



- x B 
e +i 



(43) 



This shows that Mq is a solution of the homogeneous equation with an error [the second term in (|43|l ] which vanishes 
linearly as xb —> 0, showing that Mo is accurate to 0(x B +e °). This is not very satisfactory, particularly when A is 
small and t§ — ► — 1. The pole term, not yet discussed, will correct this problem. The second term in (|42H gives a 
correction proportional to vp?xbIQ 2 ■ At large Q 2 it becomes 



lim r.h.s. 

Q 2 — oo 



2nd 



M (x b ,k,Q 2 ) - C 



(to 2 + k)x b + Q 2 (l - x%) 



e +i 



m 2 XB 



Q 2 {l + x B ) 



(44) 



In applications to DIS, these terms can safely be neglected. 

Next, we look at the effect of inserting the pole term (|25|l into the r.h.s. of 1)27(1 . This gives 



r.h.s. 



2-K\e 2 x 2 B f 1 da r ( x B 

pole _ Q2 

27rA e 2 xb 

~ Q 2 



Q 2 



[m 2 + K'} 2 



m 2 XB 



to 2 + xbk (to 2 + k) xb + Q 2 (l — x 2 B ) 



(45) 



where we assumed that Q 2 y> to 2 so that the delta function gives a = xb- The second term is small at large Q 2 
(except at the exceptional point xb = !)■ Combining the first term in i|45|) with the correction term in l|43|l gives 



r.h.s. 



pole+corr 



27rA e 2 x B 

Q 2 





'2ir\e 2 




~ XB 


[ Q 2 C0 _ 


T 2 

X B 



- c x B 
2ir\e 2 



m 2 + xbk 
-(4 + 1) c 



(46) 
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Setting 



27rAe 2 



Co 



(47) 



will cancel the leading error term and insure that the 
ansatz Ij41ll is accurate to 0(x 2 B +e °) and vanishes asis -> 
for all £q. We conclude that the solution (|41|l with 
the choice (14 7|) is accurate for all x B and Q 2 as long as 
Q 2 » to 2 . 



D. Summation of the series for small ij and n 

Before leaving this discussion, it is instructive to show 
explicitly how the Regge behavior l|40l) arises by summing 
the infinite number of terms in the series defined by the 
iteration of Eq. l(77|) . 

Iteration of Eq. (|27|> generates the following series: 

M(x b ,k,Q 2 ) - M po i c (x B ,n, Q 2 ) 

oo 

= Y J M {n+1 \x B ,K,Q 2 ) 

n=0 

where, taking K max — * oo, 

J*, K > [m 2 + k'} 2 

with 

M {0) (x Bl n,Q 2 ) = M polc (x B ,n,Q 2 ) 

and 



(48) 



(49) 



(50) 



16tt 2 



da / dn' . 

XB ^ OtK 



(51) 



Near n = 0, the leading log^s) terms in this series can 
be summed. In Appendix B we show that these leading 
terms are given by the simpler series 



M££L(x B ) = A / da M^ ding (x B /a) 



leading V 

With the aid of the identity 



1 / Xb rhi 1 

a V , ni \ 1 i n+l 

— !og (y) = — — 7 log ^ 
n y n+l 

the general term becomes 

c x B 



1 

J'b 



(n+l) 



(x B ) 



'leading ^, 

and summing to all orders gives 



A log 



1 

x B 



(52) 



(53) 



(54) 



c x B 



oo ^ 

^2 7T\ 



n=0 



A log 




fl\ 


Cb)_ 


= Cf)X B — 






\x B J 



— ( 55 ) 



in precise agreement with the small A and small k limit 
of PI) (where i ~-l + \). 

This calculation confirms the normalization the solu- 
tion, but it also shows explicitly that the Regge behavior 
arises when Alogxs ~ — 1> the condition that insures 
that all of the loops are of comparable size. Under these 
conditions perturbation theory breaks down and a non- 
perturbative summation of all diagrams is needed. 



IV. DIS STRUCTURE FUNCTIONS 

In this section the results of the previous section are 
applied to DIS. First we add the photon vector degrees 
of freedom, then we compute the two diagrams shown in 
Fig. HI 



A. Photon degrees of freedom 

The vector degrees of freedom of the photon will in- 
troduce additional structure into the pole term (|25|l and 
the scattering amplitude (|41|) . The pole term becomes 

J pole( X b,P,Q) 

= (2p - qr(2p - q) u M pole (x B ,-p 2 , Q 2 ) . (56) 

If both the initial and final quarks were on-shcll, this 
would be gauge invariant, but the initial quark is not 
on-shell. We impose gauge invariance in the manner sug- 
gested in Ref. [Rj , where it was found that the substitu- 
tion 



J» -» J" 



q 2 



(57) 



insured that the final state interactions vanished in the 
DIS limit. This choice therefore preserves the dominance 
of the pole term in the DIS limit, and is what is needed 
in the absence of a full theory. Applied to the product of 
two currents, it leads to the replacement 



{J» v q v .)q v (J""V)9' 1 



(J M v q v rq v >)t<t 



(58) 



This, in turn, give the substitutions 







2 - . 

<r 


pV- 


- (V 


p-qq^ 
q 2 


?Y- 














p-qq 

q 2 



(59) 

and the pole term Q56f) is transformed into the familiar 
form 



J^b,p, Q) = ^p u M pole (x B , k, Q 2 ) . 



(60) 
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where, as before, k — — p 2 . The general structure of the 
full amplitude will be expressed in term of two scalar 
functions 



J^{x B ,p,Q) 



—35- W2(XB,K, Q ) , 



(61) 



where the Wi differ from the familiar Wi of DIS by nor- 
malization and the dependence on the off-shell quark 
mass. 

When the form l|61[) is inserted into the integral that 
leads eventually to Eq. ll'-'il) or (|27|) , the integration over 
the azymuthal angle 7r, which was originally trivial, now 
requires use of the identity 



1 

2^ 



2n 



1.0. „„ p^p v 



p 2 



-k 2 zr + 



-hi + (po - fc ) 2 • (62) 



The quantities in this equation can be expressed in terms 
of the variables used in Eq. (|2"7|) . In the small x b , large 
Q 2 limit 



(1 — a) (k' — an) 



~2 2 

P =Po 



4*| 



1 



k 2 ± + (po - fc ) 2 
a 2 Q 2 



2 

(p~kf 



a 2 Q 2 
4x% 



(63) 



Hence Eq. I|27|l generalizes to a system of coupled equations for the Wi 



W! (x b ,k, Q 2 



16tH 



1 

da 



d.K' 



2x 2 B (1 - a) (k' - an) 



W2(xb,k,Q 2 ) - 4p 2 M po io(*B,K,(3 2 ) 



16tt 2 



da 



W2(x B /a, k', Q 2 ) + wi(xB/a, k', Q 2 )| 



dn 



, W2 {x B /a, k', Q 2 ) 
\m 2 + k'] 2 



(64) 



Note that the equation for W2 is identical to (|27|l except 
for the factor of 4p 2 multiplying the pole term, and that 
w\ is smaller than W2 by Q~ 2 , and can be ignored. The 
factor of Ap 2 changes the value of Co extracted from the 
cancellation described in Eq. I)46|) , but otherwise does not 
alter the argument. Incorporating these modifications 
into the solution obtained in the previous section, and 
taking the high Q 2 limit, gives 



wi{x B ,n,Q 2 ) 
W2{x Bl n,Q 2 ) 





2ite 2 5(x B - 1) 
1 

xb 



e 2 g 2 
8nm 2 



m 2 + k 



(65) 



Note that this solution is dimensionless and independent 
of Q 2 . 



B. Structure function for scalar bound states 

We now use the result (|65|) to construct the contribu- 
tions from each of the diagrams shown in Fig. ^ Figure 
0] shows the labeling of momenta in the generic convolu- 
tion diagram. In this part of the calculation we choose a 
different coordinate system, with momenta 

P = {A/ B ,O,0} 
q = {u,0,q L } 



{M B - £,p±,p 3 } , 



(66) 



where £ — \/M 2 



Pi 



j»3 as required by the condition 



that the heavy quark is on its positive energy mass-shell. 
Since the off-shell quark structure function u>2 is mani- 
festly covariant, it may be written in any coordinate sys- 
tem, and we will use this property to write it in the co- 
ordinate system defined by the choices given in Eq. II66II . 
With this choice the variable x takes on the familiar form 



Q 2 



2M B v 



(67) 



In the covariant spectator formalism|ll|. the diagram 
in Fig. U| is 



>p, q) 



where J^ v was defined in Eq. Ij61|l and 

T(p 2 ) 



(68) 



(69) 



is the relativistic wave function (by definition), with 
r(p 2 ) the covariant vertex function for the coupling of a 
bound meson with mass Mb to an on-shell heavy quark 
with mass M and an off-shell light quark with four- 
momentum p. In the scalar theory we are describing this 
can be a function of p 2 only. 
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0(pi,j*) 



(P-p) 



(72) 




FIG. 4: Diagram showing the convolution of the scattering 
amplitude W2 from Eq. Ij65^ with the covariant bound state 
vertex function, as canculated in the spectator theory. The 
on-shell heavy quark is labeled by an x in the diagram. 



Decomposing (1681) into the standard scalar functions 



J^(P, q) = Q 2 ) r v + W*(u, Q 2 ) 



(70) 



multiplying by and P^P^, and separating W\ and 
W2, gives the following expressions 

W x {v, Q 2 ) = J | Wl + X -w 2 [6 - 1] J , (71) 

with 



and we retained Wi for completeness. In the high Q 2 
limit, q L ^ v ^ Q 2 , and P • p — > (P ■ q)(p ■ q)/Q 2 and 
p 2 -> {p-q) 2 /Q 2 giving 



lim 0(pi,p3)->l. 

Q' 1 — »oo 



(73) 



Hence Wi vanishes (as expected for scalar constituents), 
and, as Q 2 — > 00, 

fW 2 ^, Q 2 ) - 2M s x | Mg| * V) W2 (74) 



The integral 174|l is best expressed in terms of the in- 
dependent variables p\ and p~ . Introducing the dimen- 
sionless variable z = P-/Mb, so that 



P = 



Xb 



1 - z 
r>2 



s + Q 2 z 
and introducing the constant 



(p 2 ± + zM 2 - z(l - z)M 2 B ) 
— (at large Q 2 ) , 



S = m 2 M B ^— 
An 



(75) 



(76) 



Wi = Wi(x B , -p ,Q ) 



gives 



vW 2 (v, Q 2 ) 



xS / — f 

171 



dz 
1- z 



5( X --l 



x 



m — p 



(77) 



Note that 1/W2 goes like x 1 l ° at small x, as it is expected 
on general grounds [1-4] . 

It is instructive and amusing to evaluate vW% for a 
simple model, and look at the comparative behavor of the 
two terms. A detailed study would require the solution 
of the bound state spectator equation for the spectator 
wave function. While this is not difficult, it is also not 
necessary in order to obtain a qualitative understanding 
of the physics. To this end we introduce a very simple 
ansatz for the wave function 



V(p 2 ) 
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m — p z 



M 2 (z)+p\ 



(78) 



with 



M 2 (z) = (1 - z)m 2 + zM 2 - z(l - z)M\ 



(79) 



In the following example we take j3 = 1 for simplicity. 

We now evaluate each of the terms in (|77|l . The first 
term is the pole contribution, shown in Fig.^a). Imme- 
diate evaluation of the z integral gives 



pole 



S'f) 



1-x 



dp 2 



*V(* = *))i (80) 



and evaluating the p± integral gives 

S Q x 2 (1-x) 



vWi 



pole 



1 — x + x Mfi 



o(l - x)M B 



(81) 



BO 
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This function is shown (as a dashed line) in Fig. [S] for 
the choices Mq = M/ra — 2, Mbo = Ms/m = 2, and 
5*0 = 10. It gives the familar shape anticipated in the 
drawing shown in Fig. ^ it is zero at the x = and 
x = 1 boundaries. 

The new result of this paper is the inclusion of the 
small x Regge term with a normalization consistent with 
the valance contribution. For the model H78[) this term 
becomes 



vW 2 



A 



with 



GO, A) 



Regge 



S 



dz , 



GO, A) 



m 2 (l — z) 
M 2 (z) 



(82) 



(83) 



where the power l§ was given in Eq. (|37|) . [Note that 
Eq. (l8*2*|) agrees with the general result of Eq. .] Since 
1 — ^0 < 2, the Regge term is always larger than the pole 
term near x = 0, but the importance of this term depends 
on the size of A. If Iq > (A > 2) the Regge contribution 
to the quark distribution f(x) = vW 2 /x is singular as 
x —>■ 0, and the structure function has the shape shown in 
Fig. n But, even for small values of A, Fig.[S]shows that 
the Regge term enhances the distribution significantly. 



V. CONCLUSIONS 

This paper applies the Regge concepts originally de- 
veloped in Ref. Q to DIS. We find that an approximate 
solution can be obtained that can explain, qualitatively 
at least, the enhancement of the structure function at 
small x. The physics that produces this effect is the in- 
elastic production of any number of massless gluons when 
x is small, and s is corresponding large. 

We emphasize that this calculation shows that, even 
when the strong coupling is small, so that perturbation 



1.5 



1.0 



/(*) 



0.5 







/ /// " 




/ Ir 




11/ 




Ir 





0.0 



0.2 



0.6 



0.8 



1.0 



FIG. 5: The structure function f(x) = vW%/x as a function 
of x for the model discussed in the text. The solid lines are 
the sum of the pole and the Regge contributions for values of 
A = 0.2, 0.5, 1.0 and 2.0, with the curves increasing in value as 
A increases. The dashed line is the pole contribution, which 
is independent of A in our simple approximation. 



theory should apply over most of the phase space, there 
are regions when an an infinite number of diagrams must 
be summed to get the correct result. The need to do 
this for the description of bound states is familar to most 
physicists. Here we have discussed it is also necessary 
in the small x region when the production of massless 
gluons must be summed to all orders to give the correct 
result. 

The scalar theory examined in this paper gives an un- 
realistically small result £(0) ~ — 1 + A for the t = 
intercept of the Regge trajectory. A realistic treatment 
based on QCD, including spin and confinement, should 
give an intercept of £(0) = 1/2, corresponding to the p 
trajectory (the largest trajectory possible for a qq lad- 
der). While this result could be mimicked in this model 
by choosing A = 3.75, this is too large to be realistic, and 
the failure is due to the physics missing from the scalar 
theory. 

The gluon distribution, which corresponds to £(Q) ~ 1, 
is expected to arise from single and multiple pomeron 
exchanges mechanisms not included in the model 
discussed here. The calculation of ladder exchanges of 
massless gluons of spin one is known to give an inter- 
cept around one for small coupling constant, in agree- 
ment with the small x behavior of structure functions. 

There are many interesting issues and unsolved prob- 
lems yet to be dealt with. The very notion of Regge tra- 
jectories ai(t) in QCD is understood to some extent only 
for positive t p|, where it has been shown that the trajec- 
tories are indeed linear and the intercepts aAO) = £i(0) 
have been explicitly calculated theoretically ■ At the 
same time, the behavior of cti(t) for t < in QCD (and 
other field theory models) is still unknown, though very 
interesting from both experimental and theoretical point 
of view. 

As trajectories for t < are understood, it will be 
of interest to apply these ideas to other processes, such 
as exclusive DIS and Virtual Compton Scattering, where 
t 7^ and Q 2 ^ Q' 2 . Such applications will yield mod- 
els for the interesting Generalized Parton Distributions 
(GPDs). Another point of interest is the interconnection 
of Regge behavior and the DGLAP evolution of struc- 
ture functions. All of this will require an explicit un- 
derstanding of the behavior of the Regge trajectories at 
both positive and negative values of t. To obtain lin- 
ear Regge trajectories with realistic intercepts, and to 
describe structure functions quanitatively in the whole 
x region, one needs to modify the ladder kernel of the 
equations used in this paper to include the nonpertur- 
bative effects arising from confinement. Study of these 
fundamental and difficult problems is planned for future 
work. 
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may be converted to a partial derivative equation through 
the following steps. 

First, after taking y- derivative from both sides in 
Ea. (|A3|l . we obtain 



APPENDIX A 



In this appendix we give an alternative demonstration 
that study of Eq. 1291) near x = will not produce solu- 
tions that have not already been considered. 

Rewrite Eq. (HJJ 



(Al) 



where use the fact that the upper limit of the y' integra- 
tion can be taken to infinity with errors that vanish as 
Q 2 -> oo. 

To simplify the form of the integral in this equation, 
let us define a new function <j)(y, x) trough 



r{y,x) = (f){y,x){ y + l) 2 x l + 1 . 
The equation for <f> 



(A2) 



<j>(y,x)(y + l) 2 = A 



da 



a 

x ya 



dy^(y'A) 



d(</>(y,x)(y + l) 
dy 



-A / da<j) [ya, — j . (A4) 



Changing the integration variable from a to z = — , we 
write Ea. (|A4p as 



d 



1 j:r 



where the notations 



(Hy^ x )(y + !) 2 ) = - Xx J dzcj) (xyz, , 

(A5) 



$ = 4>{y, x)(y + l) 2 = r(y, x)x~ [l+1) , $ y = — $ 

ay 



are used. 



For partial x and y - derivatives of $ y Ea. (|A5|l . we 
(A3) have 

I 



y,x 



1/x 



— § y = -$ v - [ ) ^x(j){y,x) - Xx j dz yz 



l/x 

d f 
'I',,.,, = —'!>,, A.r / d: ,r; 



i=xyz 



i=xyz 



The integral in the last term in Eq . Q A6(l is the same as in Ea. (|A7fl . therefore Eo. (|A6(l may be re- written as 



<& y ,x = -® y + -<t>{y,x) + -3> 

X X X 



y,y> 



or, finally, 



y%, v + * 



A 



" (y + i) 2 
r 



$ - X<&y. X = 0. 



(A6) 
(A7) 

(A8) 
(A9) 



The form of the & y . x ~ term in Eq. l|A9(l suggests a 
direct way to separate the variables. A simple observa- 
tion is that for a function of the form F(y,x) = f(y)x n 
with arbitrary n, the last term in Eo. (|A9|) with mixed 



derivatives becomes a term with only y - derivative: 

- xF v , x = -nF y . (A10) 
Therefore, we will look for a solution of Ea. l|A9(l in the 
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form 



®(y,x) 



-(i+i) 



(All) 



with an additional condition that it should be also con- 
sistent with the original integral equation IjAlfl . Here 



with I to be determined. Note that values of n here 
are not necessary integer but are to be found from the 
consistency with the resuling equations for r n (y). For the 
homogeneous equation i|A9|l , coefficients c n are arbitrary. 
Each r n (y) is a solution of the equation 



x 



o. 



(A12) 



L n = I + 2 - n. 



The general solution is 



(A13) 



r n (y) = d (y + l) A F(a, b, c; y+l) + C 2 (y + l) s F(a', b', c': y + 1), (AH) 

with 

a = A, a' = B, 

b = L n - 1 + A, V = L n -l + B, 

c = 2A, d = 2B, (A15) 

A= i(l + VTT4A) = 1 + 2 B= VI+4A) =1-A = -(I + 1). 

This presentation has the form most suitable for the analysis of the consistency with the additional conditions. 

The Hypergeometric series in each of the terms of (|A12|) becomes a finite sum, if b (&') is a negative integer or zero; 
and the series does not converge for positive y otherwize. For consistency with the integral equation, it is necessary 
that all r n (y) vanish as y — > oo. This leads to the following conditions: a) Ci=0, and b) b' should be a non-positive 
integer satisfying —b' < —B. 

Combining all expressions, we obtain 

r n (y) = T—^ F(-(l + l),-n,-2(l + iy,y + l), n < I + 1 = ~ (vT+4A - l) . (A16) 

For n — 0, this reproduces the solution Ea. H39|) . For sufficiently small values of A (A < 2), it is the only solution. For 
2 < A < 6, there two different solutions with n — and n = 1, and so on. We have recovered the solutions already 
discussed in Section 3. , 



APPENDIX B 

In this Appendix it is shown that Eq. (|49(l can be ap- 
proximated by the simpler Eq. (|52l) . This equation says 
that the leading xb dependence of (n + l)th order term 
is obtained from the nth order term by treating its k 
dependence as a constant, fixed by its value at n = 0. 

It turns out that the best way to show this is to ex- 
amine the structure of the first two terms in the series. 
Dropping the unnecessary Q 2 dependence (and noting 
that the pole term does not depend on k) , the first order 
term is 



(m 2 + ft') 2 

Jxb \ m + aK J a ^ a 



C x B - 



C x b , 



(Bl) 



m + xbh 

where the last is exact if k — 0. Calculation of the sec- 
ond order term reveals the pattern that will allow us to 
calculate and sum the leading small xb behavior of all 
the terms. Inserting the exact result for gives 



xb, k) 



XCqXb 



1 da 



(B2) 



dn' 



(m 2 



c l\2 



However, the k' integral is dominated by the region near 
its lower limit n' ~ aK, and in this region the factor 
arising from the k dependence of M^' is unity (up to 
terms of order xg) 



\ + 0{x B ). (B3) 



rn- 



xbk 
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Hence the leading order result for M^> is simply 
MV( XB , K ) * XC 0XB J £ £ Mj^p , 



(B4) 



precisely what would have been obtained if had 
been evaluated at k = from the start, and its k depen- 
dence ignored. Furthermore, the structure of Eq. (|B4|> 
shows that the leading x B dependence of is ob- 

tained by setting n = 0. When xb —* 0, the a integral is 
dominated by values of a near xb, and hence 



m 



m 



dn'- 



(m + xbk) 



K (m 2 + k') 2 
1 + 0(x B ). (B5) 



This confirms that the leading order contribution to 
is obtained from 

mSL^M =X[ da M^{x B /a) (B6) 



in full agreement with Eq. I|52l) . This can be confirmed 
order by order, proving l|52l) . 



We conclude by confirming by direct calculation some 
of our previous arguments. First, the exact result for the 
k' integral in (|B2|) is 



Ib2 



dn'- 



(m 2 



-M2 



X B 



log 



1 



a — xb 
m 2 (a — xb) 
a(m 2 + x b k) 



m z + n a — xb 
This is regular, even when a — > xb, and hence 

r t m2 



(B7) 



ok 



(B8) 



which is the result we obtain using the approximation 
(0. 

Next, we check the approximation used to obtain M^> 
by evaluating the integral in (|B4|I exactly: 



1 da 



m 



a m 2 + an 



= log — lo 

x B 



m + k 
m 2 + igK 



(B9) 



When xb is small, only the log(l/;EB) term survives, 
which is equal to the exact result at K = 0. This confirms 
the correctness of the approximation (|B5|) . 
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